Cortical neurons obtained from patient-derived iPSCs with GNAO1 p.G203R variant show altered differentiation and functional properties.

Pathogenic variants in the GNAO1 gene, encoding the alpha subunit of an inhibitory heterotrimeric guanine nucleotide-binding protein (Go) highly expressed in the mammalian brain, have been linked to encephalopathy characterized by different combinations of neurological symptoms, including developmental delay, hypotonia, epilepsy and hyperkinetic movement disorder with life-threatening paroxysmal exacerbations. Currently, there are only symptomatic treatments, and little is known about the pathophysiology of GNAO1-related disorders. Here, we report the characterization of a new in vitro model system based on patient-derived induced pluripotent stem cells (hiPSCs) carrying the recurrent p.G203R amino acid substitution in Gαo, and a CRISPR-Cas9-genetically corrected isogenic control line. RNA-Seq analysis highlighted aberrant cell fate commitment in neuronal progenitor cells carrying the p.G203R pathogenic variant. Upon differentiation into cortical neurons, patients' cells showed reduced expression of early neural genes and increased expression of astrocyte markers, as well as premature and defective differentiation processes leading to aberrant formation of neuronal rosettes. Of note, comparable defects in gene expression and in the morphology of neural rosettes were observed in hiPSCs from an unrelated individual harboring the same GNAO1 variant. Functional characterization showed lower basal intracellular free calcium concentration ([Ca2+]i), reduced frequency of spontaneous activity, and a smaller response to several neurotransmitters in 40- and 50-days differentiated p.G203R neurons compared to control cells. These findings suggest that the GNAO1 pathogenic variant causes a neurodevelopmental phenotype characterized by aberrant differentiation of both neuronal and glial populations leading to a significant alteration of neuronal communication and signal transduction.

HuD (ELAVL4) gain-of-function impairs neuromuscular junctions and induces apoptosis in in vitro and in vivo models of amyotrophic lateral sclerosis.

Early defects at the neuromuscular junction (NMJ) are among the first hallmarks of the progressive neurodegenerative disease amyotrophic lateral sclerosis (ALS). According to the "dying back" hypothesis, disruption of the NMJ not only precedes, but is also a trigger for the subsequent degeneration of the motoneuron in both sporadic and familial ALS, including ALS caused by the severe FUS pathogenic variant P525L. However, the mechanisms linking genetic and environmental factors to NMJ defects remain elusive. By taking advantage of co-cultures of motoneurons and skeletal muscle derived from human induced pluripotent stem cells (iPSCs), we show that the neural RNA binding protein HuD (ELAVL4) may underlie NMJ defects and apoptosis in FUS-ALS. HuD overexpression in motoneurons phenocopies the severe FUSP525L mutation, while HuD knockdown in FUSP525L co-cultures produces phenotypic rescue. We validated these findings in vivo in a Drosophila FUS-ALS model. Neuronal-restricted overexpression of the HuD-related gene, elav, produces per se a motor phenotype, while neuronal-restricted elav knockdown significantly rescues motor dysfunction caused by FUS. Finally, we show that HuD levels increase upon oxidative stress in human motoneurons and in sporadic ALS patients with an oxidative stress signature. On these bases, we propose HuD as an important player downstream of FUS mutation in familial ALS, with potential implications for sporadic ALS related to oxidative stress.

Prediction of protein-RNA interactions from single-cell transcriptomic data.

Proteins are crucial in regulating every aspect of RNA life, yet understanding their interactions with coding and noncoding RNAs remains limited. Experimental studies are typically restricted to a small number of cell lines and a limited set of RNA-binding proteins (RBPs). Although computational methods based on physico-chemical principles can predict protein-RNA interactions accurately, they often lack the ability to consider cell-type-specific gene expression and the broader context of gene regulatory networks (GRNs). Here, we assess the performance of several GRN inference algorithms in predicting protein-RNA interactions from single-cell transcriptomic data, and propose a pipeline, called scRAPID (single-cell transcriptomic-based RnA Protein Interaction Detection), that integrates these methods with the catRAPID algorithm, which can identify direct physical interactions between RBPs and RNA molecules. Our approach demonstrates that RBP-RNA interactions can be predicted from single-cell transcriptomic data, with performances comparable or superior to those achieved for the well-established task of inferring transcription factor-target interactions. The incorporation of catRAPID significantly enhances the accuracy of identifying interactions, particularly with long noncoding RNAs, and enables the identification of hub RBPs and RNAs. Additionally, we show that interactions between RBPs can be detected based on their inferred RNA targets. The software is freely available at https://github.com/tartaglialabIIT/scRAPID.

CircAFF1 Is a Circular RNA with a Role in Alveolar Rhabdomyosarcoma Cell Migration.

Circular RNAs (circRNAs), covalently closed RNAs that originate from back-splicing events, participate in the control of several processes, including those that occur in the development of pathological conditions such as cancer. Hereby, we describe circAFF1, a circular RNA overexpressed in alveolar rhabdomyosarcoma. Using RH4 and RH30 cell lines, a classical cell line models for alveolar rhabdomyosarcoma, we demonstrated that circAFF1 is a cytoplasmatic circRNA and its depletion impacts cell homeostasis favouring cell migration through the downregulation of genes involved in cell adhesion pathways. The presented data underline the importance of this circular RNA as a new partial suppressor of the alveolar rhabdomyosarcoma tumour progression and as a putative future therapeutic target.

FUS Alters circRNA Metabolism in Human Motor Neurons Carrying the ALS-Linked P525L Mutation.

Deregulation of RNA metabolism has emerged as one of the key events leading to the degeneration of motor neurons (MNs) in Amyotrophic Lateral Sclerosis (ALS) disease. Indeed, mutations on RNA-binding proteins (RBPs) or on proteins involved in aspects of RNA metabolism account for the majority of familiar forms of ALS. In particular, the impact of the ALS-linked mutations of the RBP FUS on many aspects of RNA-related processes has been vastly investigated. FUS plays a pivotal role in splicing regulation and its mutations severely alter the exon composition of transcripts coding for proteins involved in neurogenesis, axon guidance, and synaptic activity. In this study, by using in vitro-derived human MNs, we investigate the effect of the P525L FUS mutation on non-canonical splicing events that leads to the formation of circular RNAs (circRNAs). We observed altered levels of circRNAs in FUSP525L MNs and a preferential binding of the mutant protein to introns flanking downregulated circRNAs and containing inverted Alu repeats. For a subset of circRNAs, FUSP525L also impacts their nuclear/cytoplasmic partitioning, confirming its involvement in different processes of RNA metabolism. Finally, we assess the potential of cytoplasmic circRNAs to act as miRNA sponges, with possible implications in ALS pathogenesis.

The PRALINE database: protein and Rna humAn singLe nucleotIde variaNts in condEnsates.

SUMMARY: Biological condensates are membraneless organelles with different material properties. Proteins and RNAs are the main components, but most of their interactions are still unknown. Here, we introduce PRALINE, a database for the interrogation of proteins and RNAs contained in stress granules, processing bodies and other assemblies including droplets and amyloids. PRALINE provides information about the predicted and experimentally validated protein-protein, protein-RNA and RNA-RNA interactions. For proteins, it reports the liquid-liquid phase separation and liquid-solid phase separation propensities. For RNAs, it provides information on predicted secondary structure content. PRALINE shows detailed information on human single-nucleotide variants, their clinical significance and presence in protein and RNA binding sites, and how they can affect condensates' physical properties. AVAILABILITY AND IMPLEMENTATION: PRALINE is freely accessible on the web at http://praline.tartaglialab.com.

TGS1 impacts snRNA 3'-end processing, ameliorates survival motor neuron-dependent neurological phenotypes in vivo and prevents neurodegeneration.

Trimethylguanosine synthase 1 (TGS1) is a highly conserved enzyme that converts the 5'-monomethylguanosine cap of small nuclear RNAs (snRNAs) to a trimethylguanosine cap. Here, we show that loss of TGS1 in Caenorhabditis elegans, Drosophila melanogaster and Danio rerio results in neurological phenotypes similar to those caused by survival motor neuron (SMN) deficiency. Importantly, expression of human TGS1 ameliorates the SMN-dependent neurological phenotypes in both flies and worms, revealing that TGS1 can partly counteract the effects of SMN deficiency. TGS1 loss in HeLa cells leads to the accumulation of immature U2 and U4atac snRNAs with long 3' tails that are often uridylated. snRNAs with defective 3' terminations also accumulate in Drosophila Tgs1 mutants. Consistent with defective snRNA maturation, TGS1 and SMN mutant cells also exhibit partially overlapping transcriptome alterations that include aberrantly spliced and readthrough transcripts. Together, these results identify a neuroprotective function for TGS1 and reinforce the view that defective snRNA maturation affects neuronal viability and function.

A multifunctional locus controls motor neuron differentiation through short and long noncoding RNAs.

The transition from dividing progenitors to postmitotic motor neurons (MNs) is orchestrated by a series of events, which are mainly studied at the transcriptional level by analyzing the activity of specific programming transcription factors. Here, we identify a post-transcriptional role of a MN-specific transcriptional unit (MN2) harboring a lncRNA (lncMN2-203) and two miRNAs (miR-325-3p and miR-384-5p) in this transition. Through the use of in vitro mESC differentiation and single-cell sequencing of CRISPR/Cas9 mutants, we demonstrate that lncMN2-203 affects MN differentiation by sponging miR-466i-5p and upregulating its targets, including several factors involved in neuronal differentiation and function. In parallel, miR-325-3p and miR-384-5p, co-transcribed with lncMN2-203, act by repressing proliferation-related factors. These findings indicate the functional relevance of the MN2 locus and exemplify additional layers of specificity regulation in MN differentiation.

MicroRNAs and other small RNAs in Aedes aegypti saliva and salivary glands following chikungunya virus infection.

Mosquito saliva facilitates blood feeding through the anti-haemostatic, anti-inflammatory and immunomodulatory properties of its proteins. However, the potential contribution of non-coding RNAs to host manipulation is still poorly understood. We analysed small RNAs from Aedes aegypti saliva and salivary glands and show here that chikungunya virus-infection triggers both the siRNA and piRNA antiviral pathways with limited effects on miRNA expression profiles. Saliva appears enriched in specific miRNA subsets and its miRNA content is well conserved among mosquitoes and ticks, clearly pointing to a non-random sorting and occurrence. Finally, we provide evidence that miRNAs from Ae. aegypti saliva may target human immune and inflammatory pathways, as indicated by prediction analysis and searching for experimentally validated targets of identical human miRNAs. Overall, we believe these observations convincingly support a scenario where both proteins and miRNAs from mosquito saliva are injected into vertebrates during blood feeding and contribute to the complex vector-host-pathogen interactions.

Lnc-SMaRT Translational Regulation of Spire1, A New Player in Muscle Differentiation.

The destiny of a messenger RNA is determined from a combination of in cis elements, like peculiar secondary structures, and in trans modulators, such as RNA binding proteins and non-coding, regulatory RNAs. RNA guanine quadruplexes belong to the first group: these strong secondary structures have been characterized in many mRNAs, and their stabilization or unwinding provides an additional step for the fine tuning of mRNA stability and translation. On the other hand, many cytoplasmic long non-coding RNAs intervene in post-transcriptional regulation, frequently by direct base-pairing with their mRNA targets. We have previously identified the lncRNA SMaRT as a key modulator of the correct timing of murine skeletal muscle differentiation; when expressed, lnc-SMaRT interacts with a G-quadruplex-containing region of Mlx-γ mRNA, therefore inhibiting its translation by counteracting the DHX36 helicase activity. The "smart" mode of action of lnc-SMaRT led us to speculate whether this molecular mechanism could be extended to other targets and conserved in other species. Here, we show that the molecular complex composed by lnc-SMaRT and DHX36 also includes other mRNAs. We prove that lnc-SMaRT is able to repress Spire1 translation through base-pairing with its G-quadruplex-forming sequence, and that Spire1 modulation participates to the regulation of proper skeletal muscle differentiation. Moreover, we demonstrate that the interaction between DHX36 and lnc-SMaRT is indirect and mediated by the mRNAs present in the complex. Finally, we suggest an extendibility of the molecular mechanism of lnc-SMaRT from the mouse model to humans, identifying potential functional analogues.

Circular RNA ZNF609/CKAP5 mRNA interaction regulates microtubule dynamics and tumorigenicity.

Circular RNAs (circRNAs) are widely expressed in eukaryotes and are regulated in many biological processes. Although several studies indicate their activity as microRNA (miRNA) and protein sponges, little is known about their ability to directly control mRNA homeostasis. We show that the widely expressed circZNF609 directly interacts with several mRNAs and increases their stability and/or translation by favoring the recruitment of the RNA-binding protein ELAVL1. Particularly, the interaction with CKAP5 mRNA, which interestingly overlaps the back-splicing junction, enhances CKAP5 translation, regulating microtubule function in cancer cells and sustaining cell-cycle progression. Finally, we show that circZNF609 downregulation increases the sensitivity of several cancer cell lines to different microtubule-targeting chemotherapeutic drugs and that locked nucleic acid (LNA) protectors against the pairing region on circZNF609 phenocopy such effects. These data set an example of how the small effects tuned by circZNF609/CKAP5 mRNA interaction might have a potent output in tumor growth and drug response.

CircVAMP3: A circRNA with a Role in Alveolar Rhabdomyosarcoma Cell Cycle Progression.

Circular RNAs (circRNAs), a class of covalently closed RNAs formed by a back-splicing reaction, have been involved in the regulation of diverse oncogenic processes. In this article we describe circVAMP3, a novel circular RNA overexpressed in RH4, a representative cell line of alveolar rhabdomyosarcoma. We demonstrated that circVAMP3 has a differential m6A pattern opposed to its linear counterpart, suggesting that the two isoforms can be differently regulated by such RNA modification. Moreover, we show how circVAMP3 depletion in alveolar rhabdomyosarcoma cells can impair cell cycle progression, through the alteration of the AKT-related pathways, pointing to this non-coding RNA as a novel regulator of the alveolar rhabdomyosarcoma progression and as a putative future therapeutic target.

catRAPID omics v2.0: going deeper and wider in the prediction of protein-RNA interactions.

Prediction of protein-RNA interactions is important to understand post-transcriptional events taking place in the cell. Here we introduce catRAPID omics v2.0, an update of our web server dedicated to the computation of protein-RNA interaction propensities at the transcriptome- and RNA-binding proteome-level in 8 model organisms. The server accepts multiple input protein or RNA sequences and computes their catRAPID interaction scores on updated precompiled libraries. Additionally, it is now possible to predict the interactions between a custom protein set and a custom RNA set. Considerable effort has been put into the generation of a new database of RNA-binding motifs that are searched within the predicted RNA targets of proteins. In this update, the sequence fragmentation scheme of the catRAPID fragment module has been included, which allows the server to handle long linear RNAs and to analyse circular RNAs. For the top-scoring protein-RNA pairs, the web server shows the predicted binding sites in both protein and RNA sequences and reports whether the predicted interactions are conserved in orthologous protein-RNA pairs. The catRAPID omics v2.0 web server is a powerful tool for the characterization and classification of RNA-protein interactions and is freely available at http://service.tartaglialab.com/page/catrapid_omics2_group along with documentation and tutorial.

A longitudinal study defined circulating microRNAs as reliable biomarkers for disease prognosis and progression in ALS human patients.

Amyotrophic lateral sclerosis (ALS) is a fatal neurodegenerative disease associated with motor neuron degeneration, muscle atrophy and paralysis. To date, multiple panels of biomarkers have been described in ALS patients and murine models. Nevertheless, none of them has sufficient specificity and thus the molecular signature for ALS prognosis and progression remains to be elucidated. Here we overcome this limitation through a longitudinal study, analyzing serum levels of circulating miRNAs, stable molecules that are recently used as promising biomarkers for many types of human disorders, in ALS patients during the progression of the pathology. We performed next-generation sequencing (NGS) analysis and absolute RT quantification of serum samples of ALS patients and healthy controls. The expression levels of five selected miRNAs were quantitatively analyzed during disease progression in each patient and we demonstrated that high levels of miR-206, miR-133a and miR-151a-5p can predict a slower clinical decline of patient functionality. In particular, we found that miR-206 and miR-151a-5p serum levels were significantly up-regulated at the mild stage of ALS pathology, to decrease in the following moderate and severe stages, whereas the expression levels of miR-133a and miR-199a-5p remained low throughout the course of the disease, showing a diagnostic significance in moderate and severe stages for miR-133a and in mild and terminal ones for miR-199a-5p. Moreover, we found that miR-423-3p and 151a-5p were significantly downregulated respectively in mild and terminal stages of the disease. These data suggest that these miRNAs represent potential prognostic markers for ALS disease.

Epigenetic regulation of Wnt7b expression by the cis-acting long noncoding RNA Lnc-Rewind in muscle stem cells.

Skeletal muscle possesses an outstanding capacity to regenerate upon injury due to the adult muscle stem cell (MuSC) activity. This ability requires the proper balance between MuSC expansion and differentiation, which is critical for muscle homeostasis and contributes, if deregulated, to muscle diseases. Here, we functionally characterize a novel chromatin-associated long noncoding RNA (lncRNA), Lnc-Rewind, which is expressed in murine MuSCs and conserved in human. We find that, in mouse, Lnc-Rewind acts as an epigenetic regulator of MuSC proliferation and expansion by influencing the expression of skeletal muscle genes and several components of the WNT (Wingless-INT) signalling pathway. Among them, we identified the nearby Wnt7b gene as a direct Lnc-Rewind target. We show that Lnc-Rewind interacts with the G9a histone lysine methyltransferase and mediates the in cis repression of Wnt7b by H3K9me2 deposition. Overall, these findings provide novel insights into the epigenetic regulation of adult muscle stem cells fate by lncRNAs.

Intronic Determinants Coordinate Charme lncRNA Nuclear Activity through the Interaction with MATR3 and PTBP1.

Chromatin architect of muscle expression (Charme) is a muscle-restricted long noncoding RNA (lncRNA) that plays an important role in myogenesis. Earlier evidence indicates that the nuclear Charme isoform, named pCharme, acts on the chromatin by assisting the formation of chromatin domains where myogenic transcription occurs. By combining RNA antisense purification (RAP) with mass spectrometry and loss-of-function analyses, we have now identified the proteins that assist these chromatin activities. These proteins-which include a sub-set of splicing regulators, principally PTBP1 and the multifunctional RNA/DNA binding protein MATR3-bind to sequences located within the alternatively spliced intron-1 to form nuclear aggregates. Consistent with the functional importance of pCharme interactome in vivo, a targeted deletion of the intron-1 by a CRISPR-Cas9 approach in mouse causes the release of pCharme from the chromatin and results in cardiac defects similar to what was observed upon knockout of the full-length transcript.

Zooming in on protein-RNA interactions: a multi-level workflow to identify interaction partners.

Interactions between proteins and RNA are at the base of numerous cellular regulatory and functional phenomena. The investigation of the biological relevance of non-coding RNAs has led to the identification of numerous novel RNA-binding proteins (RBPs). However, defining the RNA sequences and structures that are selectively recognised by an RBP remains challenging, since these interactions can be transient and highly dynamic, and may be mediated by unstructured regions in the protein, as in the case of many non-canonical RBPs. Numerous experimental and computational methodologies have been developed to predict, identify and verify the binding between a given RBP and potential RNA partners, but navigating across the vast ocean of data can be frustrating and misleading. In this mini-review, we propose a workflow for the identification of the RNA binding partners of putative, newly identified RBPs. The large pool of potential binders selected by in-cell experiments can be enriched by in silico tools such as catRAPID, which is able to predict the RNA sequences more likely to interact with specific RBP regions with high accuracy. The RNA candidates with the highest potential can then be analysed in vitro to determine the binding strength and to precisely identify the binding sites. The results thus obtained can furthermore validate the computational predictions, offering an all-round solution to the issue of finding the most likely RNA binding partners for a newly identified potential RBP.

Proteomics analysis of FUS mutant human motoneurons reveals altered regulation of cytoskeleton and other ALS-linked proteins via 3'UTR binding.

Increasing evidence suggests that in Amyotrophic Lateral Sclerosis (ALS) mutated RNA binding proteins acquire aberrant functions, leading to altered RNA metabolism with significant impact on encoded protein levels. Here, by taking advantage of a human induced pluripotent stem cell-based model, we aimed to gain insights on the impact of ALS mutant FUS on the motoneuron proteome. Label-free proteomics analysis by mass-spectrometry revealed upregulation of proteins involved in catabolic processes and oxidation-reduction, and downregulation of cytoskeletal proteins and factors directing neuron projection. Mechanistically, proteome alteration does not correlate with transcriptome changes. Rather, we observed a strong correlation with selective binding of mutant FUS to target mRNAs in their 3'UTR. Novel validated targets, selectively bound by mutant FUS, include genes previously involved in familial or sporadic ALS, such as VCP, and regulators of membrane trafficking and cytoskeleton remodeling, such as ASAP1. These findings unveil a novel mechanism by which mutant FUS might intersect other pathogenic pathways in ALS patients' motoneurons.

Trans-generational epigenetic regulation associated with the amelioration of Duchenne Muscular Dystrophy.

Exon skipping is an effective strategy for the treatment of many Duchenne Muscular Dystrophy (DMD) mutations. Natural exon skipping observed in several DMD cases can help in identifying novel therapeutic tools. Here, we show a DMD study case where the lack of a splicing factor (Celf2a), which results in exon skipping and dystrophin rescue, is due to a maternally inherited trans-generational epigenetic silencing. We found that the study case and his mother express a repressive long non-coding RNA, DUXAP8, whose presence correlates with silencing of the Celf2a coding region. We also demonstrate that DUXAP8 expression is lost upon cell reprogramming and that, upon induction of iPSCs into myoblasts, Celf2a expression is recovered leading to the loss of exon skipping and loss of dystrophin synthesis. Finally, CRISPR/Cas9 inactivation of the splicing factor Celf2a was proven to ameliorate the pathological state in other DMD backgrounds establishing Celf2a ablation or inactivation as a novel therapeutic approach for the treatment of Duchenne Muscular Dystrophy.

Modulation of circRNA Metabolism by m6A Modification.

N6-methyladenosine (m6A) is an RNA modification well-known for its contribution to different processes controlling RNA metabolism, including splicing, stability, and translation of mRNA. Conversely, the role of m6A on the biogenesis and function of circular RNAs (circRNAs) has yet to be addressed. circRNAs belong to a class of covalently closed transcripts produced via a back-splicing reaction whereby a downstream 5' splice donor site fuses to an upstream 3' splice acceptor site. Starting from circ-ZNF609 as a study case, we discover that specific m6As control its accumulation and that METTL3 and YTHDC1 are required to direct the back-splicing reaction. This feature is shared with other circRNAs because we find a significant direct correlation among METTL3 requirement, YTHDC1 binding, and the ability of m6A exons to undergo back-splicing. Finally, because circ-ZNF609 displays the ability to be translated, we show that m6A modifications, through recognition by YTHDF3 and eIF4G2, modulate its translation.